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Recently developed Generalized Cosmic Chaplygin gas (GCCG) is studied as an unified model of 
dark matter and dark energy. To explain the recent accelerating phase, the Universe is assumed to 
have a mixture of radiation and GCCG. The mixture is considered for without or with interaction. 
Solutions are obtained for various choices of the parameters and trajectories in the plane of the 
statefinder parameters and presented graphically. For particular choice of interaction parameter, we 
have shown the role of statefinder parameters in various cases for the evolution of the Universe. 

PACS numbers: 



I. INTRODUCTION 

Recent observations of type la Supernovae indicate that the expansion of the Universe is accelerating 
[1-5] and lead to the search for a new type of matter which violates the strong energy condition, i.e., 
p + 3p < 0. The matter responsible for this condition to be satisfied at some stage of evolution of 
the universe is referred to as dark energy [6-8]. Several candidates to present dark energy have been 
suggested with observations: the cosmological constant [7, 9], quintessence [10, 11], phantom [12, 13], 
braneworld models [14], pure Chaplygin gas model [15], generalized Chaplygin gas (GCG) model [16, 
17], modified Chaplygin gas (MCG) model [19, 20]. In the GCG and MCG approach dark energy and 
dark matter can be unified by using an exotic equation of state (EOS). Interesting feature of MCG (or 
GCG) EOS is that it shows radiation era (or dust era) in the past while a ACDM model in the future. 

In 2003, P. F. Gonzdlez-Diaz [21] have introduced the generalized cosmic Chaplygin gas (GCCG) model 
in such a way that the resulting models can be made stable and free from unphysical behaviours even 
when the vacuum fluid satisfies the phantom energy condition. The EOS of this model is 

P=-p- a [C + ( P 1+a -C)~ w ] (1) 

where C — — 1 with A a constant which can take on both positive and negative values and 

— / < w < 0, I being a positive definite constant which can take on values larger than unity. The EOS 
reduces to that of current Chaplygin unified models for dark matter and dark energy in the limit w — > 
and satisfies the conditions: (i) it becomes a de Sitter fluid at late time and when w = — 1, (ii) it 
reduces to p — wp in the limit that the Chaplygin parameter A — ► 0, (iii) it also reduces to the EOS of 
current Chaplygin unified dark matter models at high energy density and (iv) the evolution of density 
perturbations derived from the chosen EOS becomes free from the pathological behaviour of the matter 
power spectrum for physically reasonable values of the involved parameters at late time. This EOS 
shows dust era in the past and ACDM in the future. 



Since models trying to provide a description of the cosmic acceleration are proliferating, there exists the 
problem of discriminating between the various contenders. To this aim Sahni et al [22] proposed a pair 
of parameters {r, s}, called statefinder parameters. In fact trajectories in the {r, s} plane corresponding 
to different cosmological models demonstrate qualitatively different behaviour. The above statefinder 
diagnostic pair has the following form: 
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where H (= -) and q (= — fy) are the Hubble parameter and the deceleration parameter respectively. 
The new feature of the statefindcr is that it involves the third derivative of the cosmological radius. These 
parameters are dimensionless and allow us to characterize the properties of dark energy. Trajectories 
in the {r, s} plane corresponding to different cosmological models, for example ACDM model diagrams 
correspond to the fixed point s = 0, r = 1. 



In this paper, we consider the Universe is filled with the mixture of radiation and GCCG in section 
II. We perform a statefinder diagnostic to this model without and with interaction in different cases 
in sections III and IV respectively. From statefinder parameters we have shown graphically that 
the universe starts from radiation era instead of dust era. Different phases of the evolution of the 
universe have been shown graphically. With interaction case, the model goes from radiation to 
ACDM era only and without interaction case the model goes from radiation to ACDM and further 
from ACDM to phantom era and then back to ACDM. The paper ends with a short discussion in section V. 



II. MIXTURE OF GCCG AND RADIATION 

The metric of a spatially flat isotropic and homogeneous Universe in FRW model is, 

ds 2 = dt 2 - a 2 (t) [dr 2 + r 2 (d6 2 + sin 2 9d<p 2 )] (3) 
where a(t) is the scale factor. 

The Einstein field equations are (choosing 8itG = c = 1) 

3— = ptot 4 



and 



6^ = -(ptot + 3ptot) (5) 



The energy conservation equation {T" = 0) is 

Ptot + S ~(Ptot + Ptot) = (6) 
where, p to t and p to t are the total energy density and the pressure of the Universe, given by, 

Ptot = P + Pr (7) 

and 

Ptot =P + Pr (8) 

with p and p are respectively the energy density and pressure due to the GCCG satisfying the EOS (1) 
and p r and p r are the energy density and the pressure corresponding to the radiation fluid with EOS, 

Pr = IPr (9) 

where 7=4- 

Since GCCG can explain the evolution of the Universe starting from dust era to ACDM, considering 
the mixture of GCCG with radiation would make it possible to explain the evolution of the Universe 
from radiation to ACDM. 
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III. WITHOUT INTERACTION 



In this case GCCG and the radiation fluid are conserved separately. Conservation equation (6) yields, 



and 



p + 3-(p + p) = 



Pr + Z-{p r +Vr 



From equations (1), (9), (10), (11) we have 

p= C + ( 1 + 

and 



B 



3(l+a)(l+tu) 



-3(1+7) 



p r = p a 

For the two component fluids, equation (2) takes the following forms: 
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r = 1 + 



2(p + Pr ) 



dp, , dp r , ' 

d-p ip+p) + dp^S Pr+Pr \ 



and 



(P+Pr) 

Also the deceleration parameter q has the form: 



dp . , dp r . . 



1 3 (p+Pr 



1 aH 2 2 2 Vp + Pr 
Now substituting u = y = — , equation (14) and (15) can be written as, 



r = 1 



and 



2(1 + y) 



, C (u-C)- w \ (aC a, 
1 ^ '- < + -(u- 



s = 



\ u u 

2(r-l)(l + y) 



(u-C)- w +w(l + a){u-C)- 



1 



C (u-C)- 
7W — — — '— 



(10) 



(11) 



(12) 



(13) 



(14) 



(15) 



(16) 



+ 7(1 + l)v 
(17) 



(18) 



Normalizing the parameters we have shown the graphical representation of the {r, s} parameters in 
figure 1. From the figure we have seen that the universe starts from radiation era (r = 3, s > 0) via dust 
stage (2.3 < r < 2.4, s ±oo) to ACDM (r = 1, s = 0) model for choices of C = 1, B = 1, a = 1, w = 
-2,po = l- 



IV. WITH INTERACTION 

We consider the GCCG interacting with radiation fluid through an energy exchange between them. 
The equations of motion can be written as, 

p + 3% + p) = -3H5 (19) 



w - 
-2.5; 

-5 ; 
-7.5; 




r 

Fig.l 



Fig. 1 shows the variation of s against r (eqs. (17) and (18)) for C = 1, B = 1, a = 1, w = —2, po = 



and 



pV + 3-(p r +p r ) = 3ff<5 



where 8 is a coupling function. 
Let us choose, 



5 = e 



(p 1+a -C)- w 



Now equation (19) together with equation (1) gives, 

p = \c + (l - e + Bo 3-(i+°0(i-H«0) * 

Also equations (9), (20) and (22) give 

p r = Po a- 3(1+1 »+3 £fl - 3 ( 1 ^/ 

with 



/ = 



1 



dx (x 1+w +e-l) (i+™)(i+°> 



B 



and 



35(1 + 0)7 (C + x) — 

\-e + Ba-^ 1 + w ^ 1+a '> 



x = 



From (22), we see that if e = 0, i.e., 5 = 0, then the expression (22) reduces to the expression (12 
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Now for the two component interacting fluids with equations of motion (19) and (20), the {r, s} 
parameters read: 



r = 1 + 



2{p + Pr ) 



(26) 



and 



2(r- l)(p + p r 



Also the deceleration parameter q has the form: 



(27) 



1 f, , , } P + Pr 

q= o 1 + 3— Z — 

2 \ p + p r 

Now substituting u — p 1+a , y = equation (14) and (15) can be written as, 



r = 1 + 



2(1 + 2/) 



|( i+ H) +7 { (i+7)y ^ 



(28) 



(29) 



and 



where, 



2(r-l)(l + y) 

9(?+7?y) 



(30) 



C+ (l-6+Ba 3 -( 1+t> )( 



i_ 



j, = £» a -3d+7) + 3 1 a -3(l+7)j 
P P 



- = --{c + ^-cr™) 

p u 1 



and 



<5 _ (u-C)- 
p u 



£ = ^ + % + «,(! + «)(« - C) — 1 



Now we find the exact solution for the {r, s} parameters for the following particular choices of w. 



(i) If — ( 1+ ^uT+al — 1 = 0, i-e., w = \+ a " ? equation (23) can be written as 



P, = Poa- 3(1+7) -|«- 3(1+7) P (31) 

as I = — jg(c + x) 1 + a . Normalizing the parameters, the corresponding statefinder parameters are given 
in figure 2. From the figure we have seen that the universe starts from radiation era (r — 3, s > 0) via 
dust stage (2.5 < r < 2.6, s — ► ±oo) to ACDM (r = 1, s = 0) model and further from ACDM to phan- 
tom era (r < 1, s > 0) and then back to ACDM for choices of C = 1, B = 1, a = 1, w = — §, Po = 1, e = f ■ 



G 





(ii) If — ( 1+ ^,)"(i+ a ) — 1 = 1, i.e., w = 2(1+0)° ' ec l ua ti on (23) can be written as 

Pr = po a" 3 ^ - !il-ll a -3(i+7) £a ' 3(1+7) ^x 2+ - 2 Fi[2 + — , 3 + -£] (32) 

F F B 2 B 2 (l + a){2 + w)C— 1 '1 + a' C l 1 ' 

Normalizing the parameters, the corresponding statcfinder parameters are given in figure 3. 
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Fig.4 

Fig. 4 shows the variation of s against r (case (iii) )forC = l, B = l,a = 1, po = l,e = ^. 



From the figure we have seen that the universe starts from radiation era (r = 3, s > 0) via 
dust stage (1.9 < r < 2, s — > ±oo) to ACDM (r = 1, s = 0) model for choices of C = 1, 
B = l,a = l,w = -f ,po = l,e = §■ 



(iii) If w = —2, equation (23) can be written as 



a -3(i+ 7 ) 



p r = p a a 3(1+7) -t— — t r 1+2a _ 7 — AppellFi 



(l + 2a-7) 2 ,i^ (7 i+o 



1 + 2a — 7 a a — 7 2 + 3a — 7 a; 

x — £e 



(1 + a) '1 + a' (1 + a)' (1 + a) ' C : 

Normalizing the parameters, the corresponding statefinder parameters are given in figure 4. From the 
figure we have seen that the universe starts from radiation era (r = 3, s > 0) via dust stage (2.3 < r < 
2.5, s ±00) to ACDM (r = 1, s = 0) model and further from ACDM to phantom era (r < 1, s > 0) 
and then back to ACDM for choices of C = 1, B = 1, a = 1, po = 1, e = 5. 



V. DISCUSSION 



In this work, we have considered the matter in our Universe as a mixture of the GCCG and radiation 
as GCCG can explain the evolution of the Universe from dust era to ACDM. These gases are taken both 
as non-interacting and interacting mixture. In the first case we have considered a non-interacting model 
and plotted the {r, s} parameters. As expected this model represents the evolution of the Universe from 
radiation era to ACDM with a discontinuity at dust. In the second case the interaction term is chosen in 
a very typical form to solve the corresponding conservation equations analytically. Also the statefinder 
parameters arc evaluated for various choices of parameters and the trajectories in the {r, s} plane are 
plotted to characterize different phases of the Universe. These trajectories show discontinuity at same 
r in the neighbourhood of r = 2 and have peculiar behaviour around r = 1. The {r, s} curves have two 
branches on two sides of the asymptote. The branch on the right hand side of the asymptote corresponds 
to decelerating phase before (or up to) dust era, while the left hand side branch has a transition from 
decelerating phase upto ACDM era. Some peculiarity has been shown in figures 2 and 4 around r = 1. 
In these two cases, the model goes further from ACDM to phantom era and then back to ACDM. 
Moreover, in figure 4, there is further transition from ACDM to decelerating phase and then then again 
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back to ACDM. Thus we can conclude that the present model describes a number of transitions from 
decelerating to accelerating phase and vice- versa. 
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